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Abstract 

We study the sensitivity of the dilepton production cross section to higher twist terms. 
A method for calculating the resummed cross section is introduced in moment space, where 
the leading higher twist terms can be easily parametrized and calculated. We find that a 
1/Q power correction may significantly affect the cross section at fixed-target energies. 



1 Introduction 



Large perturbative corrections are common in hadron-hadron scattering cross sections in 
QCD [|I|. In this paper, we will discuss corrections associated with soft gluons that produce 
large threshold effects. In the case of dilepton production cross sections, such large correc- 
tions arise when most of the partonic momenta is carried off by the lepton pair. We will 
exhibit a class of related nonperturbative corrections that can also have a significant impact 
on this cross section. 

Several methods for resumming threshold corrections have been developed 0. Recently, 
resummed cross sections for dilepton production at fixed-target energies have been calculated 
P, ^] in principal value resummation (PVR) In this scheme, the large corrections 

exponentiate into leading and next-to-leading exponents, Ei and E^l, respectively. These 
exponents contain the effects of the one-(i?L) and two-loop (Enl) anomalous dimensions and 
running coupling, and are each defined by a principal value integration. Moreover, because 
singularities associated with the infrared (IR) behavior of the QCD running coupling are 
avoided, no explicit IR cutoffs need to be imposed in principal value resummation. An 
analysis of perturbative resummation ||^, 0] shows the presence of nonperturbative power 
suppressed (higher twist) contributions beginning at 0{1/Q). 

It was found in [^] that for protons on fixed-targets, the resummed hard part (in PVR) 
when combined with parton distributions from global fits led to cross sections that overesti- 
mate the data by roughly 30 to 50%. It was also observed in IQ that the discrepancy could 
be due to higher twist effects. The aim of this paper is to explore this possibility. 

In order to analyze the sensitivity of principal value resummation to higher twist effects, 
we introduce a technique for calculating cross sections which exploits the exponentiation of 
large corrections in moment-space. We shall show below that to a good approximation, the 
partonic flux may be approximated by a polynomial in z = Q'^/s, with Q the pair mass and 
s the partonic center-of-mass invariant mass squared. In this method, the resummed cross 
section becomes a simple and finite sum of exponentials with coefficients that are functions of 
energy alone, and most importantly, in which the higher twist effects are easily parametrized 
and calculated. 

The paper is organized as follows: in section 2, we describe in detail how we calculate 
the resummed cross section in moment space. Here, we will focus only on mass-rapidity 
distributions at zero rapidity, although our technique can be applied to other distributions 
as well. In section 3, we present numerical results and compare with experiment. We 
will find that the resulting resummed cross section in moment space, combined with terms 
proportional to M/Q, is sensitive to M of the order of hundreds of MeV. Indeed, we find 
good agreement with data when M ~ — 1 GeV, for which the higher twist contribution 
effectively cancels the effects of perturbative resummation. This, in effect, confirms that the 
cross section may be quite sensitive to higher twist, at least for fixed-target energies. Finally, 
in section 4, we give a summary and conclusions. Some numerical details are discussed in 
an appendix. 



2 



2 The cross section in moment- space 

Schematically, the dilepton-production reaction is given by 

hi{pi) + h2{p2)-^n{Qn + x, (1) 

where is the lepton-pair mass squared. We will use the standard notation s = {pi +^2)^ 
and r = Q'^/s. 

At fixed-target energies, where the Z boson may be neglected, principal value resumma- 
tion provides the following expression (in the DIS scheme), for the differential cross section 
in dilepton production at zero rapidity, 



dQ'^dy 



y=0 



^B{Q')T.^f dzufj{z,a) ff ' ' \ (2) 

J- JT Z 



with 



osiQ-) ^ ^g, (3) 



uj{z,a) = A{a)(^S{l - z) 



^ . sin[7rPi(^, «)]r(l + Pi(^, a)) 



7r(l - z) 



(4) 



T]f{Tlz) = Ff/,^{^z,Q)Fj/,^{^z,Q). (5) 

The sum in eq. (§) is over all active quark and anti-quark fiavors, is the (electromagnetic) 
fine structure constant, a is proportional to the QCD running coupling (see eq. ([T0|) below) 
and Ff/h{x,Q) denotes the parton distributions. The function A{a), which involves the 
resummation of large, ^-independent Sudakov terms, is given explicitly in appendix A, while 
we define ^ 

Pi(x,a) = —E(x,a). (6) 
dx 

The function E{z,a), which exponentiates the large corrections, is, as mentioned above, 
given as a sum of leading and nonleading contributions. In moment-space, it takes the form 

E{n,a) = E{n,a)L + E{n,a)NL- (7) 
The leading exponent is given explicitly by 

E{n,a)L = a{g['h^-gi'h2), (8) 

with 

Ii{n,t) = 2I{n,t/2) - I{n,t), 

I{n,t) = t|^dc(^^^^)ln(l + (l/t)ln(l-C)), 

- //<nr^) i + (i/t)in(i-c) ' 
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where 

a = «s(g')/7r, t = l/(«62), (10) 

with 

b2 = {llCA-2nf)/12. (11) 

It will be sufficient to consider only the leading exponent E{n, a) l in the following discussion. 
For completeness, however, we have reproduced the explicit form of the next-to-leading 
exponent E{n,a)NL, as well as the constants gi^\g2^^ (eq. (||)) in appendix A. 

We are now ready to discuss the higher twist corrections to the resummed cross section 
(P(j/dQ'^dy\y=o (eq. (||)). To do so, it will be necessary to go to moment- or n-space, where 
the resummed hard part of eq. M) takes the simple form 



1 

71-1, 



<j 

= A(a)e^("'"\ (12) 



with E{n, a) given by eqs. (0-|D and A{a) by eq. ( |3lD of appendix A. 

The leading higher twist term implied by the resummed exponent E[n, a) is of the form 

n-. (13) 

This can be derived in the following way |0. Consider the leading exponent E{n, which 
depends on Ji and I2 (eq. (^). By a change of variables, ln(l — C) ^' 1 and using 

(1 _ - 1 = V ilzZ^^-, (14) 

where (a)^ = V{a + m)/V{a) is the Pochammer symbol, I{n,t) in eq. (|]) can be rewritten 
as 

I(n, t)=ty ^ -^V / c/x'e"^" ln(l + x'/t), (15) 

where V denotes a principal value integration. After several more simple steps, one obtains 

J(n,t) = -f:^%^(mt), (16) 

where the function S is proportional to the exponential integral function: 

S{x) = xe-^P r dy—. (17) 
J-oc y 

Similarly, it can easily be shown that 



°° (I -n) 
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The leading higher twist term of E{n, a) comes from the 2agi^^I{n, t/2) term that appears 
in eq. (H). Taking only the leading (m = 1) term in the sum in eq. ([T6|) , we get 

2ag\'hin, t/2) ~ [n-jV dy-, (19) 

where we have used the 1-loop form of a, so that t (eq. (p!0D ) is simply, t = InQ^/A^. Since 
the pole of the integrand is at ?/ = 0, independent of Q, the basic structure of the leading 
power correction to E{n, a) isl/Q. This argument is specific to principal value resummation, 
but the result is quite general, as discussed in where it was also shown that this power 
correction is actually induced by the leading infrared renormalon ||^. 

We now digress on the issue of including higher twist terms in momentum- or z- space. 
With the correspondence 

n ^ (20) 
the leading power correction can be parametrized by making the substitution 

E{z,a)^E{z,a)--^^, (21) 

I - zQ 

for some constant c, in eq. (H). We then find that the direct calculation of the resummed 
cross section given by eq. (g) is hopeless, because an essential singularity in ojfj{z, a) (eq. (^)) 
will be generated aX z = 1. This highlights one of the advantages of working in ra-space, 
for, as we shall now see, the higher twist term (eq. (|TB|)) is in a form that makes it easy to 
calculate its effects in the resummed cross section. 

To write the cross section in terms of moments, we first need to express the sum involving 
the parton luminosity polynomial in z: 

Y: e}j^^fir/z) ^ Cn{T)z- + A^(r, z), (22) 
/ "=i 

where AjF denotes the error in the polynomial approximation. Some details on how rimax 
is determined, as well as how the interpolating polynomial is obtained, are described in 
appendix B. 

Using eq. ( p2|) in eq. (|]), we obtain 

A2 rimax „i „i /\T(t 7^ 

—— ^asiQ^) E Cn{T) / dzz--^ujfj{z,a) + aB{Q^) / dzufj{z,a) ' \ (23) 

"V "1/ y=0 j^—i Jr Jr Z 

Combining eqs. ([T^) and (|23|) , and applying several straighforward manipulations, we 
find 



f^max 



- ^BiQ')A{a) Y Cn(r)e^('^'°) + CT^ + CT^. (24) 

The explicit forms of the correction terms CTi and CT2 will be given shortly. To study the 
sensitivity of eq. ( p4D to higher twist, we insert a term proportional to nA/Q in the exponent, 
so that now the n-space resummed cross section is given by 



dQ^dy 



c^aBmAia) Y C„(r)exp 
y=o n=l 



E{n, a) — 



cnA 



+ CTi + CT2, (25) 
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where A is the QCD scale, c is a constant, and the correction terms are 



^1 ^0 7r(l - z) 

withQ cujj(z, a) given by eq. @, and 

^^(l:£).^43Z^_"f:'c„M.-. (27, 



We will use eq. (^) to calculate resummed cross sections in the next section. Note that in 
this form, the cross section is simply a finite sum of exponentials weighted by the coefficients 
of the interpolating polynomial, up to corrections (CTi and CT2) which are expected to be 
very small, since AjF(r, z) is small. 



3 Numerical results 

For our numerical calculations, we chose to work with the kinematics of the E605 experiment 



T0[] , which studied proton beam interactions with a copper target at ^/s =38.8 GeV. We 
have used the CTEQ2D distributions with A at 4(5) fiavors equal to 0.235(0.155) GeV. We 
have also included all finite 1-loop contributions (see ||^) in all the resummed cross sections 
that we are now going to present. In the discussions to follow, we will simply denote the cross 
section (fa/dQ^dy\y=o by a. Hence, a{nmax, cA) will refer to the cross sections computed in 
n-space from eq. (^5]) as described above, for a given Umax and cA, and cr^ the cross section 
computed by the usual method as an integral over z (eq. (0)) with no explicit higher twist, 
at the same value of Q. We also point out that the cross section defined by eq. ( pSj) is 
dominated by the sum over n, with the correction terms giving a contribution of not more 
than 1.6%. 

To determine rimax, "we used the method described in [jll|. The basic idea is that by adding 



a z"'^^ term to a degree-n polynomial fit to a particular set of data (the parton luminosity 
(eqs. (^, 13))), the deviation of the interpolation from the data is decreased. However, for 
some degree rimax, the addition of higher powers of z no longer decreases the deviations by 
a statistically significant amount. We have found that for the data points with which we 
compare, Umax ranges from 13 to 16 only. We also mention that, although the choice of Zmin, 
which defines the range of validity of the interpolation (see appendix B), somewhat affects 
the value of Umax, we have verified that the n-space resummed cross section is insensitive to 
That is, a{nmax,cA) ~ a{n'^^^,cA), where 

'^maxi'^'max) obtained when ZminiZmin) 

used in the interpolation. 

We plot in fig. 1 the resummed cross sections calculated in moment-space and scaled 
relative to the corresponding z-space cross sections taken from the resummed curve (solid 
line in fig. 2). The points denoted by circles, triangles and squares correspond to n-space 

"'^Note that for brevity, we have represented the 1/(1 — 2:) dependence of the functions in CT2 as simply z. 
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cross sections calculated using eq. ( pSj) with cA equal to 0, 0.5 and 1 GeV, respectively [recall 
that we take c = for cr^] . 

We see that the n-space resummed cross section (T{nmax,0) is different from the ^-space 
cross section az- The two differ, however, by no more than 12%, with cr{nmax,0) initially 
smaller but eventually becoming larger than az- This difference already illustrates the range 
of sensitivity to higher twist effects that are implicit in any resummation prescription. We 
can also infer from fig. 1 that for cA = 1(0.5) GeV, the higher twist effects are about 27 to 
23% (14 to 11%) of the resummed cross section cr{nmax,0), decreasing with r. 

In fig. 2, we plot the data of experiment E605 and reproduce the curves of fig. 9d of 
0] (dashed curve=2-loop and solid line=resummed cross section). For comparison, we also 
include the results obtained for a^rimax, IGeV) (dotted curve). We see that, with this value 
of cA, the n-space resummed cross section gives quite a good fit to the data for this range of 
r. In fact, the a^rimax, IGeV) curve is almost identical to the 1-loop curve in fig. 9d of [Q. 

4 Discussions and conclusions 

In this paper, we have derived an expression in moment-space for the cross section in principal 
value resummation. In this form, the cross section becomes a finite sum of exponentials, and 
higher twist terms can be included in an easy, straightforward manner. This is in contrast 
with the resummation formula written in momentum- or z-space, where it is not possible to 
treat the higher twist in an analogous fashion. 

Direct comparison with experiment shows that a. 1/Q correction, with a coefficient of 
order 1 GeV, in this variant form of principal value resummation, agrees with experiment, 
while the purely perturbative resummed cross sections in |Q overestimate the data. It seems 
that such a higher twist term is large enough to cancel the effects of resummation, bringing 
down the resummed cross section to the size of the 1-loop cross section. We also note that 
this size of the higher twist is comparable to those found in e'^e~ event shape variables ||12|| . 

We also mention that for the particular distribution we are concerned with, (Pa/ dQdy\y=o, 
it was found in that already at 0{as), the non-singular terms in the hard part account 
for about 20% of the total cross section. Although we have included all such 1-loop finite 
contributions to our resummed cross sections, there is still some uncertainty due to similar 
non-singular contributions at higher orders. However, since we are mostly concerned with 
the structure and effects of the leading higher twist term, and not the precise calculation of 
cross sections, our conclusions are not affected by this uncertainty. 

It should be noted that we have not really explained why the higher twist coefficient 
should be as big as 1 GeV. It could be that it is actually smaller, in which case, the resummed 
cross section will still overestimate the data but less than the amount found in 0. The 
remaining discrepancy could then be explained by the other factors alluded to in [Q: the 
importance of higher order effects when combining deeply inelastic scattering (DIS) and 
hadron-hadron data in global fits, the theoretical uncertainty in A{a), and the effects of 
finite higher order terms in the hard part. Probably, the full cross section comes from a 
combination of all such effects. Here, however, we have shown that higher twist does play a 
role. 
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Appendix A The functions E{n,a)NL and A{a) 

Here, we write down the exphcit expressions for the functions E{n, a)^^ and A{a) which 
were described in section 2. The next-to-leading exponent E{n,a)NL is given by B 



(1) 



,(2) 



,(2) 



(28) 



with 



Jl 
J2 

K2 



-(a&3/&2 
dC 



dC 



1 



P 

n-1 



dC 



1-C 

^n— 1 



C dy ln(l + (l/t)ln[(l-C)(l-y)]) 
1-y (l + (l/t)ln[(l-C)(l-t/)])2^ 
l\ln(l + (l/t) In(l-C)) 



1-C y (l + (l/t)ln(l-C))2 
C dy 1 



1-C 

^n-l _ I 



l-y{l 
1 



[i/t)H{i- Oil -yW 



C ;(i + (i/t)in(i-c)) 



2' 



(29) 



where a and t are defined in eq. (|10]). All integrals are evaluated over a principal value 
contour as defined in H, |^ . The various constants are given by0 



.(1) 



(2) 
92 



gi ' = 2Cf, g2 



(1) 



3^ „{2) _ ^ /67 7r2\ 5ny' 

2<-'F, 5'l - (^18 ~ ~ ~9~. 



- 12C(3)) + CACpCii,' - If + 6C(3)) + (-fvr^ + f ) 
63 = (34C1 - (lOC^ + 6Ci.)n^)/48, 



(30) 



where n/ is the number of fiavors, and Ca-, Cp are color factors. The parameter 62 is defined 
by eq. (Q). 

We also point out here that in the calculation of -Etvl for this paper, the term proportional 
to the integral K2 was not included. This term, strictly speaking, yields logarithms beyond 
nonleading order, as defined in ^. This treatment differs from the calculations in 
where the same term was kept. However, we have checked that the contribution from the 
above K2 term is no more than 1% of the total cross section. In particular, the resummed 
curve in fig. 2, which excludes the contributions from the K2 term, is nearly identical to the 
solid curve of fig. 9d in [||. 

The explicit form of A(a), to the same order, is 0] 



A{a) 



1 + 2CF0i + ba^ exp -Cp 



a 



^-3 
3 



Inr + 7rG(2) 



sin^ 



-a 



^(1) 1 ^(2) ^ 

+ ^-{7rb29a - In r) + In r 



26^ a 



261 



(31) 



. The relation between the g'f^ 



^The parameter g^^'' can be found from Table 1 of jl^, where i?^^^ = 
and JS*^ is described in 11 . The extra factor of 1/4 in eq. ( ^o|) is due to our expansion in a = ctg/'^ rather 
than as/2TT as in fl^ . 
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with 



1 + (vrfeatt)^ 



16 3^^ ^ 6 , 



CaCf 



2545 
432 



216 6^^ ^ 



2Cf , = C^Ci. 



--C(2) 
18 { 



11C(2) - ^C(3) 



(32) 



and 
b 



^2 



-23 



720 



-IT 



35 
96' 



-TT 



+nfCF 



rC(3) 



19 
72 



7 



+ CaCf 



144 216^ 12^^ ' 1440 



TT 



-TT 



L3'' ' 72 54 
where C(s) is the Riemann Zeta function. 



(33) 



Appendix B Interpolating Polynomials 

The interpolating polynomial that we require is defined by eq. (|22|) , repeated below: 

f n=l 

For a particular value of we generate a set of values of the left-hand side of the above 
equation in a certain range of z, Zmin < z < 1. The parameter Zmin is selected such that 
J^]f{r/z) ~ for z < Zmin- The typical values of Zmm that we used vary from 0.075 to 0.2. 
We then use the Fortran 77 IMSL package on each set of values to obtain the interpolating 
polynomial, or equivalently, the coefficients C„(r). The polynomials that we obtain fit the 
luminosity data very well, with percentage differences of the order of 10~^ or better. 

In principle, given a set of data points, one can always fit a polynomial of degree Umax-i 
up to Umax = number of points. However, as described in section 3, a method exists for 
determining a non-trivial rimax- One expects that the polynomial can fit the luminosity data 
very well in the region z > Zmin- Outside this fit region, care must be taken that the selected 
polynomial does not behave wildly. To be specific, the lowest value that z (or Zmin) can 
take is r, where the luminosity vanishes. It must then be checked that for z < Zmin, the 
polynomial fits are free of large oscillations. We have verified that the fits we have obtained 
are indeed well-behaved in this region of z. 
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FIGURE CAPTIONS 



Figure 1. Fractional deviations of the resummed cross section in n-space, 

o^inmax.cK) (eq. (p5D), from the purely perturbative cross section, (eq. (0)). 
cA = (circles); cA = 0.5 GeV (triangles); cA = 1 GeV (squares). 

Figure 2. Comparison with E605 data. 

Solid=resummed a^; Dashed=2-loop; Dotted=resummed airimaxicK = 1 GeV). 
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